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1. INTRODUCTION 

Industrial applications required high mechanical accuracy such as computer-numerical-control 
machines, robotics, electric vehicles utilizing permanent magnet synchronous motors (PMSMs) because of 
their high power, low moment of inertia, high start-up torque [1]. Similarly to induction motor drive (IMD), 
field-oriented-control (FOC) and direct torque control (DTC) strategies [2], [3] are used in high-performance 
torque and flux controls of PMSM drive [4], [5]. The DTC has a simpler structure than FOC, therefore its 
simulation and real implementation are easier than FOC. This strategy also gives high robustness and fast 
response of torque [1], [5]. 

Conventional DTC methods bring high torque and current ripples because they utilize look-up tables 
for desired voltages to control the flux and torque [3], [5]. Besides that, they make switching frequency of 
power converters not constant, even though it was fixed. Integrating pulse-width modulation (PWM) into 
DTC makes switching frequency constant, and provides improved performance of the torque and the current 
[6]-[17]. Popular PWM methods are sinusoidal PWM (SPWM) [6]-[8] and space vector PWM (SVPWM) 
[8]-[10], [12]-[14]. The SVPWM offers a better dc-link utilization of 15% compared to the SPWM [7]. 

In PMSM drive, measurement stator current errors affect performance of DTC strategies [18]. 
kalman filter (KF) and its variations have many applications listed in [19]. KF were employed to observe 
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speed and load [20], estimate currents and their time derivatives [21], [22], identify flux and inductances 
variation [23], estimate parameter combinations [24]. Forgetting factor was employed to compute noise 
covariance in dual adaptive KF algorithms [20]. Time derivatives of stator current were provided by KF in 
model reference adaptive system-based sensorless IMD [21]. Stator currents were estimated by extended KF 
(EKF) in sensorless PMSM drive [22]. Least squares method was utilized together with EKF to obtain 
PMSM parameters [23]. Two KF strategies were employed to estimate online any identifiable electrical 
parameter combinations [24]. Future observations and error minimization of variables were estimated by 
EKF for controller design [25]. In the paper, estimations and identifications are omitted, instead of that KF is 
employed to provide approximation of stator currents distorted by Gaussian measurement noises. 


2. MATHEMATICAL MODEL OF DTC-PMSM DRIVE 
Figure 1 shows control structure of DTC PMSM drive with PWM and KF. In (1)-(4) describe 
mathematical model of PMSM in [a, f] coordinate system: 


dWsa ` 
Usa = EE + iggRs (1) 
dps ; 
Usp = — + isg Rs (2) 
Wsa = Lsisa + Py cos 0, (3) 
Psp =L; isg + Wy sin 6, (4) 


Where usa, Usg are stator voltage vector components; Ya, Ys are stator flux vector components; isa, isg are 
stator current vector components; 9, is rotor position; Rs is stator resistance; Ls is stator inductance; wy is 
magnetic flux of the PM. In (5) calculates motor torque T, from stator flux and stator current components, 
and mechanical relation of motor torque, load torque Tz, and mechanical speed @,, is expressed by (6): 


Te = 1.5p (isg sa = isaWsp) (5) 
dwm _ (Te-T1) 
dt Jm (6) 


Where p is number of pole pairs, and Jm is motor inertia. 
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Figure 1. Control structure of DTC PMSM drive with PWM and KF 


Signal computation block estimates stator flux components, stator flux vector magnitude, orienting 
angle y, motor torque utilizing reconstructed stator voltages and Kalman filtered stator currents, according to 


(7-01): 
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Bsa = J (Usa — Rsîsa)dt (7) 
Psp = S (usg — Rsîsg)dt (8) 
Ds = |2 + 2, (9) 
? = arcsin ($£) (10) 
T, = 1.5p(Psaisg — Ispisa) (11) 


Where symbol ^ denotes estimated quantities. KF block smooths stator current components distorted by 
Gaussian measurement noises [19]. Next section describes the SVPWM and the KF. 


3. COMPUTATION PROCESS OF SPACE VECTOR PWM AND KF 

In order to obtain the average value of space vector with desired Usa, usg components, PWM 
techniques are utilized to provide switch-on and switch-off durations of power electronic devices of the 
inverter. For the SVPWM technique, in one switching period ts, durations ta, tp, to that are respectively time 
intervals of using vectors Ua, Us, Uo (Vo or V7) in Figure 2(a), are computed by (12)-(16): 


U = (Ua)? + (ušp)? (12) 


a = tan (=) (13) 
te =t; {cos asin 5) — sin & cos C=) aa (14) 
tp =t, {- cosa sin [e21] + sin q cos [e] aa (15) 
to = ts — ta — tp (16) 


Where sx is sector (1 to 6). Figure 2(b) shows that SVPWM utilizes better DC link voltage Uz. than SPWM, 
so the SVPWM is implemented in this paper. Switching combination in Table | is employed to reduce loss 
due to change state of power electronic devices. 


A B SVPWM 
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Figure 2. SVPWM method [7] (a) phase vector & reference trajectory and (b) comparison with SPWM 
method on maximum linear control voltage 
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Table 1. Switching combination of SVPWM [7] 
2 3 


Sector 1 4 5 6 
Upper switches S\=tyttytto/2 S\=t,t+to/2 S\=to/2 S\=to/2 S\=ty+to/2 S\=tattytto/2 
(Si, S3, Ss) S3=fytto/2 S3=totty+to/2 S3=tattytto/2 S3=tatto/2 S3=fo/2 S3=to/2 

S5=to/2 Ss=to/2 Ss=tytto/2 Ss=tattytto/2 Ss=tottytto/2 Ss=tatto/2 
Lower switches S4=to/2 Sy=ty+to/2 S4=tattytto/2 S4=tettytto/2 S4=tatto/2 S4=to/2 
(Sa, S6, S2) So=tatto/2 So=to/2 So=to/2 So=tytto/2 So=tattytto/2 So=tattytto/2 
So=tattytto/2 So=totty+to/2 So=ty +fo/2 So=fo/2 Sa=to/2 So=tytto/2 


The KF block brings estimate of stator current isa, iss components according to (17)-(24): 


xg = Gy (17) 
Yr = Hx_ + Vk (18) 
Zk = cent (19) 
P, = FP, FT (20) 
Ze = Ve — HŠę (21) 
Ky = B,H™ (HB.HT + R)“ (22) 
Ry = Žk + Kye (23) 
P,, = (I — KH) P, (24) 


Where: x=[isa isg]’: state vector; F, H: state transition matrix, measurement matrix; v: zero-mean Gaussian 
measurement noise vector with covariance R=o°I; symbol~denotes predicted vectors. 


4. RESULT AND DISCUSSION 

Parameters of PMSM are listed in Table 2. Results are implemented in MATLAB/Simulink 
environment at @n re=150 rpm, load torque jump= 0.5Ty (see Figure 3), values of covariance o°={0.25, 0.57, 
1.0°, 2.07, 4.07} (see Figures 4-5). Frequency converter has DC-link voltage of 372 Vdc, and switching 
frequency of 20 kHz. Performance of the proposed control structure and its conventional version without KF 
is evaluated by indices including ripples in forward operation (durations 0.2 s-0.3 s and 0.4 s-0.5 s) and 
reverse operation (durations 0.7 s-0.8 s and 0.9 s-1.0 s), the integral of the absolute value of speed error 
(IAE), the integral of time multiplied by the absolute value of speed error (ITAE) [26]-[28]: 


1 
IAE = |\e,,(t)|dt (25) 
0 


1 
ITAE = {tle,,(0)|dt (26) 
0 


Figures 4-5 show motor speeds in different cases of noise covariance o°. It is easy to see that the 
larger parameter o°, the higher ripple of motor speeds for both control structures. Table 3 indicates that the 
structure with KF brings reduction of 34%-81%, 53%-75% in motor speed ripple for forward and reverse 
operations respectively compared with the structure without KF, and the larger parameter o? is, the more 
ripple tends to decrease. The ripples are only computed at the times when maximum value of absolute of 
speed error happens. Table 4 respectively show non-time-related and time-related performance indices IAE 
and ITAE. The structure with KF reduces 0.04%-64%, 0.06%-63% IAE, ITAE values than the structure 
without KF. Above evaluations show that the proposed structure owns much lower performance indices than 
the structure without KF. Reason for this is that filtered stator currents leads not only to smoother stator 
current responses, but also smaller ripples of stator fluxes and motor torque (see Figures 6-8). With large 
values 2.0°, 4.0? of o, ripples of stator currents, stator fluxes and motor torque tend to increase, and it 
requires other approaches to lower the ripples, IAE, and ITAE. 
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Table 2. PMSM parameters [4] 
Symbol Quantity Value 
ny nominal speed 3000 rpm at f=200 Hz 
Ty nominal torque 7.73 Nm 
Py nominal power 2.29 kW 
In motor inertia 0.00151 kgm? 
Uinrns induced line-to-line voltage 263 V at 3000 rpm 
Isy nominal stator current 5.6A 
Wu magnetic flux of the PM 0.1706 Wb 
P number of pole pairs 4 
R, stator resistance 0.65 Q 
L, stator inductance 7.7 mH 
T T 4 7 oo I, Fr FEFE ESEA ten =a EAT =| T 
á ! ! 
2 
1 S2 4 : 
! S i i 
E i i 
aig ta ! 3 1 bons > : P ! ! a 
200 i — i i 0 a tj 
0 01 02 03 O4 05 06 07 08 09 l 0l 02 03 04 05 06 OF 08 09 l 
Time [s] Time [s] 
(a) (b) 
Figure 3. Reference (a) motor speed and (b) load torque 
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== =! = e 
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Figure 4. Motor speeds at (a) o°=0.257, (b) o°=0.5?, (c) o°=1.0?, and (d) o°=2.0? 
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Figure 5. Motor speeds at o°=4.07 
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Table 3. Motor speed ripple [rpm] Table 4. IAE and ITAE performance indices 
Forward ‘ Forward f 
: Reverse operation 2 Reverse operation 
2 operation 2 operation 
7 Without With Without yi KF ° Without With Without With 
KF KF KF KF KF KF KF 
0.25? 0.43 0.28 0.52 0.24 0.257 5.425 5.423 2.568 2.567 
0.57 1.50 0.39 1.48 0.47 0.57 5.585 5.449 2.655 2.582 
1.07 6.03 1.26 6.01 1.45 1.07 6.343 5.595 3.057 2.662 
2.0° 22.75 5.03 21.46 6.13 2.0° 10.055 6.308 4.926 3.070 
4.07 90.35 16.97 94.55 27.61 4.07 27.792 9.739 13.579 4.953 
30 r ' 30 - ; 
< 20- z 20 
-2 10- 3 
E o 8 
3 3 
= -104 ; 5 -10 
g —lsa,without KF g 
BD -20+ —"so,with KP a -20 
30 1 i i i f i 1 Z tsoywith KF 30 i L i i i i i 
“O Ool 02 03 04 05 06 07 08 09 1 “0 Ool 02 03 04 05 06 07 08 09 1 
Time [s] Time [s] 
Figure 6. Stator currents isa at (a) o°™=1.0° and (b) o°=2.07 
0.2 02 
Jsa without KF } —Vsa,withoutKF | 
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Figure 7. Stator fluxes Wa at (a) o°= 1.07 and (b) o°=2.07 


= 
© 
i 
= 
© 


n 


Torque [Nm] 
io) 

Torque [Nm] 
© 


-54 —f,, withoutK F -5 Te withoutK E 
~~" Te withK F Te with F 
-10| et -10} eletr 
(0) 0.2 04. 0.6 0.8 1 0 0.2 0.4 | 0.6 0.8 1 
Time [s] Time [s] 


Figure 8. Torques at (a) o°=1.0? and (b) o°=2.0? 


5. CONCLUSION 

The drive structure utilizing Kalman filtration and SVPWM of DTC PMSM drive was presented in 
the paper. Simulation results were carried out at a low reference speed, half of nominal torque, and large 
values of measurement noise covariance. The proposed structure dedicated significantly smaller performance 
indices including ripples, IAE and ITAE than the conventional drive structure without KF, especially at large 
covariances. It also provides smoother responses of stator current, stator flux, and motor torque compared to 
the structure without KF. Adaptive KFs can be employed to receive high-performance filtering. The proposed 
structure can be utilized in sensor or sensorless speed-control systems of PMSM drive using intelligent 
control, adaptive control or robust control. 
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